Geographical Coarsegraining of Complex Networks 
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We perform the renormalization-group-like numerical analysis of geographically embedded com- 
plex networks on the two-dimensional square lattice. At each step of coarsegraining procedure, the 
four vertices on each 2x2 square box are merged to a single vertex, resulting in the coarsegrained 
system of the smaller sizes. Repetition of the process leads to the observation that the coarsegrain- 
ing procedure does not alter the qualitative characteristics of the original scale- free network, which 
opens the possibility of subtracting a smaller network from the original network without destroying 
the important structural properties. The implication of the result is also suggested in the context 
of the recent study of the human brain functional network. 

PACS numbers: 89.75.Hc, 89.75.Fb, 05.10.Cc 
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Study of complex networks has been one of the most 
active research areas not only in physics but also in other 
various disciplines of natural and social sciences P, Q, ■ 
In some existing networks, the computerized automatic 
data acquisition techniques make it possible to grab the 
detailed information of interconnections in networks. In 
contrast, in many biological and social networks, the 
complete network structure is hard to be defined and 
even when it is possible it requires tremendous time- 
demanding efforts. The detailed structure of the neuronal 
network of Caenorhabditis elegans , composed of about 
300 neuron cells and 14 synaptic couplings per neuron, 
has been obtained by biologists through direct observa- 
tions. In comparison to the C. elegans neural network, 
the complexity of the human brain is gigantic: It con- 
tains about 10^^ neuron cells, each of which is connected 
to lO'^-lO^ other neurons via synaptic couplings. Con- 
struction of the detailed map of all neuron connections 
in human brain is beyond imagination and will be so in 
the future. 

In the viewpoint of statistical physics, on the other 
hand, understanding the qualitative collective behavior 
of the brain, although it originates from the actual de- 
tailed interactions of neurons and abundant biochemical 
substances, may not require such detailed microscopic 
map of interneuron connections. In this regard, the re- 
cent study by Eguiluz et al. in Ref. draws much inter- 
est: Brain activity has been measured from 32 x 64 x 64 
sites (called voxels) and the intervoxel correlation has 
been used to map out the functional network of the hu- 
man brain. Although the number of voxels in Ref. U 
is more than a million, each voxel still contains O(IO^) 
neuron cells. Consequently, one can say that the brain 
functional network in Ref. has been based on heav- 
ily coarsegrained information, and thus it is not clear 
whether the observed scale-freeness of the network is a 
genuine emerging property of actual interneuron connec- 
tions or not. Very recently it has been found that a very 
simple vertex merging process results in scale-free net- 
works, regardless of the initial network structure 0. In 
the present context, this observation may suggest that if 



the network is too much coarsegrained, one cannot trust 
the resulting scale-free distribution since it may not re- 
flect the structure of the original network but is a simple 
artifact of coarsegraining. 

In this Letter, we start from model scale-free networks 
that are geographically embedded 0, and then repeat 
several steps of geographic coarsegraining. We find that 
the coarsegraining process does not change important 
properties of the original network. In particular, the 
degree exponent 7, the clustering property, the assorta- 
tive feature, and the hierarchical structure do not change 
much upon the iteration of the geographic coarsegrain- 
ing. Our result suggests that the scale-free feature of the 
human brain functional network may not be the artifact 
of the coarsegraining, and thus the increase (or the de- 
crease) of the size of voxels is expected not to change the 
main results of Ref. 0. We also suggest that one can use 
the geographic coarsegraining method presented in this 
Letter to subtract a smaller network from the original 
larger network, without destroying important structural 
properties. This can be very useful when the network is 
too big to be handled for a given computational capabil- 
ity. 

We first build the geographically embedded scale-free 
network following Ref. |y: N ^ L x L vertices are put 
on lattice points of the two-dimensional square net, and 
then the degree k of each vertex is chosen according to 
the degree distribution function p{k) cx k~'^ . A vertex v 
is selected at random and then its assigned degree ky 
is realized on the basis that the geographically closer 
vertices (within the distance proportional to y/k^) are 
connected first. As the procedure is repeated over ver- 
tices, some vertices may not fulfill their assigned degrees 
if their all possible target vertices already exhausted their 
allowed number of edges. When this happens, those ver- 
tices have degrees different from the initially assigned 
ones, and there appears the cutoff degree scale (and the 
corresponding cutoff length scale) beyond which p{k) de- 
viates from the power-law form p{k) ^ k^'' (see Ref. |^ 
for details). 

Once the network is constructed in this way, we re- 
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FIG. 1: Coarsegraining procedure. Four vertices in each 2x2 
box is merged to a single vertex. After merging, two edges 
(thin solid lines) connecting vertices in different boxes be- 
comes one edge (thick solid line) with the weight w — 2. This 
procedure is repeated for all 2 x 2 square boxes resulting in 
the coarsegrained network which is four times smaller than the 
original network. We then remove edges of smaller weights in 
order to fix the average degree. 



peat the following geographic coarsegraining, which is in 
parallel to the KadanofF block spin renormalization group 
procedure in standard statistical mechanical systems (see 
Fig.^l: Four vertices on each square box of the size 2x2 
is merged to a single vertex and accordingly the edges 
connecting intra-box vertices (dashed lines in Fig.^ are 
disregarded, but the inter-box connections are kept (thin 
solid line in Fig. We also keep track of how strong 
the edges are by assigning the weight that is sim- 
ply the number of edges connecting two merged vertices 
V and w. For example, in Fig. ^ there exist two edges 
(thin solid line) connecting the two square boxes before 
the merging, which gives rise to the weight w = 2 for the 
edge (thick solid line) connecting the two vertices after 
the merging. 

If we keep all the edges in the coarsegrained network, 
the average degree increases as the procedure is iterated, 
resulting in the fully-connected network eventually. To 
remedy this, we fix the average degree at each step of 
coarsegraining by removing weaker edges with smaller 
values of the weight. Suppose that we have to remove 
Mr edges to keep the average degree the same, and that 
there are Af„, edges of the weight w. For example, for 
Mr < Mu,=i, randomly picked Mr edges oi w = 1 is 
removed. If Mw=i < Mr < Mu,=2, all w; = 1 edges re- 
moved and Mr — M^^i edges with w = 2 are randomly 
deleted . The above procedure makes sense since in real 
situations, it is common that the coarsegraining is often 
accompanied by the change of the sensitivity of the mea- 
surement: When the system is looked at a far distance, 
we only have interest in large scale structures. 

Figure |2] shows the result of the coarsegraining. Orig- 
inal networks of the size 64 x 64 are generated following 
Ref. for 7 = 2.5,3.0,4.0,5.0 and the above explained 
coarsegraining process is iterated. The network size in 
this work is smaller than the cutoff length scale beyond 
which the network ceased to be scale free, which is also 
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FIG. 2: Degree distribution p{k) versus the degree k for 
geographically embedded scale-free networks. The original 
64 X 64 networks in (a) with the degree exponents 7 = 
2.5,3.0,4.0, and 5.0 are coarsegrained n times; (b) n — 1, 

(c) n = 2, and (d) n = 3. Clearly shown is that the degree 
exponent does not change upon the iteration of coarsegrain- 
ing. Full lines in (a)-(c) are for the power law distributions 
with the exponents 2.5 and 5.0, while in (d) only the line for 
the exponent 2.5 is shown for comparison. 7 values in (b)- 

(d) only indicate the degree exponents for the corresponding 
initial networks. 
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(b) Y=3.0, 64 X 64 
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(d) 7=3.0, 128x128 



FIG. 3: Degree distribution p{k) versus the degree k for the 
networks with the degree exponent 7 = 3. The original net- 
work sizes are (a) 32 x 32, (b) 64 x 64, (c) 96 x 96, and (d) 
128 x 128. The resulting coarsegrained network at the nth 
iteration displays the same degree exponent. 



seen in Fig. Eta) where the cutoff degree scale is absent 
(see Ref.la). One sees clearly that the coarsegraining pro- 
cess does not change the degree exponent 7. In the ter- 
minology of the renormalization group (RG) formalism, 
the scale-free network with any value of 7 is the stable 
fixed point of the RG flow. This observation implies that 
the scale-free network in Ref. 6 possesses neither the de- 
gree scales nor the length scales upto the cutoff length 
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FIG. 4: The clustering coefficient C{k) versus k for the 
128 X 128 network with 7 = 3 at the n-th iteration step of 
coarsegraining for n = (a) , (b) 1, (c) 2, and (d) 3. The 
qualitative feature remains the same upon coarsegraining pro- 
cedure. 



scale which is larger than the network size in the present 
study. In Fig. |21 as the coarsegraining process is iterated 
the scale-free region appears only for sufficiently large k 
regions. 

We then study finite-size effects in Fig.|3| which shows 
p{k) at the n-th iterations for initial networks of vari- 
ous sizes (a) 32 x 32, (b) 64 x 64, (c) 96 x 96, and (d) 
128 X 128 (all for 7 = 3). Clearly exhibited is that as 
the initial network size is increased the network still re- 
mains to be scale-free even after many steps of itera- 
tions, which then excludes the possibility that observed 
behaviors are finite-size artifacts. The scale-free degree 
distribution detected in the human brain functional net- 
work ^ does not actually imply that the neural network 
of human brain is scale-free. One reason is because the 
technique in Ref.^only measures the functionality corre- 
lation of two separate voxels, not the actual path of vox- 
els through which biochemical signal transfers. One can 
also argue that since each voxel contains large number of 
neurons [about 0(10''^)], the observed scale-free distribu- 
tions can be the artifact of the coarsegraining, consider- 
ing the recent study in Ref. 5 that scale- free distributions 
can emerge from merging. Our main results 0| in the 
present study implies that this is not the case and that 
accordingly the scale-free distribution in human brain 
functional network is expected to be the genuine prop- 
erty of the brain, not the artifact of the coarsegrained 
information. 

We next investigate other important structural proper- 
ties of networks. Many real networks including Internet, 
World Wide Web, and the actor network, are character- 
ized by the existence of the hierarchical structure d, ' 
which can be usually detected by the negative correla- 
tion between the clustering coefficient (see Ref. |3) and 
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FIG. 5: (a)-(c) Density plot of the degrees fc„ and k^ where 
V and w are two vertices connected by each edge. Brighter 
region indicates that there are more edges in that region. 
The initial network shows disassortative behavior (higher fc„ 
prefers lower k^ , and vice versa) , which remains qualitatively 
the same as the coarsegraining procedure proceeds n times: 
n = (a) (initial network), (b) 1, and (c) 2. (d) Assortativity 
coefficient r versus the number n of iterations of coarsegrain- 
ing. After the coarsegraining, the network remains to be dis- 
assortative. 



the degree For example, the Barabasi- Albert net- 
work which does not possess hierarchical structure, 
is known to have the clustering coefficient C„ of the ver- 



[i.e., C(fc) - fco, see 
has been shown to 



tex V independent of its degree k 
Ref. 0] , while Holme-Kim model 

have Cik) - fc-^ ^ , in accord with the observations of 
many real networks [8j. In Fig. ^ we plot C(fc) at the 
n-th iteration step of the coarsegraining for the initial 
network of the size 128 x 128 with 7 = 3. The geograph- 
ically embedded network in Ref. is found to be some- 
how special since C{k) is better described by C(fc) ^ k^^ 
rather than the abundantly found C(fc) ^ k~^. But this 
feature remains the same upon the coarsegraining proce- 
dure, implying that the coarsegraining does not change 
the hierarchical structure of the network. 
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We next study the assortative mixing characteris- 
tics ^3 01 of the network. For the assortative network, 
vertices with the higher degree tend to have high-degree 
neighbor vertices, while for the disassortative network, 
higher degree vertices favor to have fower degree neigh- 
bors. The degrees ky and fc^ of the two vertices v and w 
connecting each edge is measured and then the histogram 
is computed by using 20 x 20 bins in log-log scales in 
fcu, plane. The brightness of the region in Fig. [3 (a)-(c) 
is chosen in proportion to the logarithm of the height of 
the histogram in that region. Again found is that the 
coarsegraining procedure does not change the disassorta- 
tive mixing property of the network, i.e., at any iteration 
step, the high-degree vertices in the network tend to have 
low-degree neighbors. This behavior of the disassortative 
mixing can also be detected by the assortativity coeffi- 
cient r (see Ref. for the definition). If r has positive 
value, the network has assortative mixing property while 
it is disassortative otherwise. In Fig. EJd), r is shown 
to have negative values at the n = 0, 1, 2 coarsegraining 
steps. The decrease of r with n is not completely under- 
stood, although this dependence of r versus the network 
size N appears to be consistent with Ref. Il^ where r 
tends to approach zero from below as the larger disassor- 
tative network is considered. 

So far we have introduced a geographical coarsegrain- 
ing procedure and applied it to the geographically embed- 
ded scale-free network in Ref. 0. Although the network 
sizes become smaller as the coarsegraining procedure pro- 
ceeds, it has been found that several key features of the 
initial networks do not change qualitatively. In particu- 
lar, the degree exponent 7 does not change, and hierar- 
chical structure [detected by the negative correlation be- 
tween the clustering coefficient C(fc) versus degree k] and 
the disassortativity (detected by that more edges connect 
high-degree vertices to low-degree vertices than to high- 
degree vertices) are remain qualitatively the same. Our 
geographic coarsegraining procedure can be useful when 
the initial network is of a huge size since one can then 
systematically reduce the network size without destroy- 
ing important characteristics of the network. Modifica- 
tion of the present coarsegraining method to apply for the 
network which is not geographical embedded can be an 
interesting extension. The main results also suggest that 
the scale-free distribution found recently for the human 
brain functional network may not be the artifact due to 
the large voxel size (each contains O(IO^) neuron cells), 
but the genuine property of the brain. 

We finally study the two-dimensional Watts-Strogatz 
(WS) network, built similarly to Ref.|^ Vertices are put 
on the two-dimensional square lattice points and every 
vertex is connected to its nearest and next-nearest neigh- 
bor vertices. Each edge is visited once, and with the 
rewiring probability P, is rewired to a randomly chosen 
vertex. The resulting network belongs to the so-called 
exponential network since the tail in the degree distri- 




FIG. 6: Geographically embedded WS network at the 
rewiring probability P — 0.1. (a) Degree distribution p{k) 
at the iteration steps n = 0, 1, 2, and 3. (b) Clustering coef- 
ficient C(n) versus n. 



bution is exponentially small. We then iterate our geo- 
graphic coarsegraining procedure with the average degree 
kept constant at each iteration. In Fig. Ela), the initial 
network of the size 128 x 128 at the rewiring probability 
P = 0.1 is coarsegrained n times. As n becomes larger, 
the degree distribution remains to be exponential and 
tends to saturate. In Fig. (SJb), initial two-dimensional 
WS networks of various sizes L = 128, 256, 512 and 1024 
at P = 0.1 are coarsegrained and the clustering coeffi- 
cient C{n) is plotted as a function of the number n of 
iterations. As the coarsegraining process proceeds the 
clustering coefficient is shown to decrease towards zero, 
which indicates that the RG stable fixed point of the WS 
network is close to the random network of Erdos and 
Renyi. 
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